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We construct, numerically, a solution of the SU(2) Bogomolny equations corresponding to a sheet 
of BPS monopoles. It represents a domain wall between a vacuum region and a region of constant 
energy density, and it is the smoothed-out version of the planar sheet of Dirac monopoles obtained 
by linear superposition. 

PACS numbers: 11.27.+d, ll.10.Lm, 11.15.-q 



O 
O 

(N 
O 

Q 



> 
o 

(N 

o 



43 



X 



I. INTRODUCTION 

Bogomolny-Prasad-Sommerfield (BPS) monopoles 
have long been of considerable interest (for reviews, see 
P, and recently have found a new interpretation in 
the context of D-branes. This D-branc connection is 
partly responsible for a particular interest in periodic 
assemblages of monopoles. For example, there have been 
studies of monopole chains Q, where the underlying 
theory (such as the Nahm transform) is well-developed. 
However, for monopole sheets or walls, where the fields 
are doubly-periodic, much less is known. 

It has been suggested 0, Q that the (essentially 
Abelian) homogeneous solution on K 3 may be viewed as 
a monopole sheet, and there is a sense in which such an 
interpretation is meaningful. If, however, one constructs 
a sheet of Dirac monopoles, as a double series, then the 
resulting field looks rather different. The purpose of this 
note is to present numerical evidence for the existence of 
an SU(2) monopole sheet which is a smoothed-out ver- 
sion of the Dirac sheet. As we shall see, it resembles a 
domain wall between a vacuum region and a region of 
homogeneous phase with constant energy density. 



II. THE U(l) CASE: SHEET OF DIRAC 
MONOPOLES 

Let r = (x, y, z) denote the usual position vector in M 3 , 
and let P be the square lattice in the xy-plane consisting 
of points of the form r = (j, k, 0), with j, k G Z. Suppose 
we place a unit-charge Dirac monopole at each site of this 
lattice. The resulting magnetic field is, at least formally, 



B(x,y,z) 



E 



(x -j,y 



k, z)/r 3 



(1) 



where r| fc 



[x — j) 2 + (y — k) 2 + z 2 . The component 



B z is a positive series which converges absolutely, and 
the other two components B x and B y are conditionally 
convergent (except, of course, at the points of P, where 



B is singular). For fixed (x,y), we have B — > (0, 0, ±7r) 
as z — > ±oo; this was checked by numerical summation of 
the series, and can also be established by a rough analytic 
approximation. 

Similarly, we may obtain a gauge potential A by taking 
a gauge potential for each monopole and summing these; 
for example, 



A(x,y,z) 



E 



(k - y,x-j,0) 
2rjk(z + rjk) 



(2) 



which is valid in the region z > 0. The asymptotic be- 
haviour of flU) is A — > J(— y, x, 0) as z — > oc. 

In general, if we have a doubly-periodic vector field B 
which is smooth for \z\ > c, then we can define its charge 
(N-,N+) by 



N± = -!- / dx j dy B z (x,y,±c); 
i7r Jo Jo 



(3) 



by continuity, N± do not change if c is increased. If B 
is genuinely a doubly-periodic magnetic field, in other 
words if its gauge potential is doubly-periodic up to a 
gauge transformation, then iV_ and N + will both be inte- 
gers: this is the Dirac quantization condition for a mag- 
netic field on the torus T 2 . The field |T|), in view of its 
asymptotic behaviour, clearly has charge {—h, h), and so 
it does not quite satisfy this quantization condition. (Of 
course, taking an array of monopoles of charge 2 would 
have the effect of multiplying all the expressions by 2, and 
so this would give a genuinely doubly-periodic magnetic 
field.) 

Another example of a 'half-doubly-periodic' field is the 
homogeneous solution 

y, z) = (0, 0, tt), A{x, y, z) = x, 0), (4) 

which has charge is (|, h). If we add (TJJ to (JTJ) , then we 
obtain a doubly-periodic magnetic field of charge (0, 1), 
with 



B 



2tt as z 
as z 



(5) 
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This is the prototype of our monopole sheet: it separates 
a vacuum region with energy density approximately zero 
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(for z <C —1), from a region of homogeneous phase with 
energy density approximately constant (for z > 1). In 
what follows, we will show that it has a non-Abelian ver- 
sion in which the singularities are smoothed out. 



III. THE SU(2) CASE: SETUP 



Let {Aj(x, y, z), &(x,y, z)} denote an SU(2) Yang- 



Assuming for simplicity that N + > > N_, and using 
Stokes's theorem, this leads to the inequality 

E L > 8n 2 L(Nl +Nl), (8) 

with equality if and only if the Bogomolny equations 



(9) 



Mills-Higgs field, and define 



dp 



are satisfied. 

There is an exact solution of (|9|) representing a field of 



Bj ■= ^s jk i(d k Ai-diA k + [A k ,Ai}), as usual. We impose charge (1,1) — in other words, with N + = 1 



-N- 



the global conditions: 

(a) Aj and $ are 2x2 anti-Hermitian trace-free ma- 
trices, and are smooth on M 3 ; 

(b) Aj and $ are periodic in both x and y (actually 
periodic, not merely up to a gauge transformation), 
with unit periods. 

In addition, we need boundary conditions as z — > ±oo, 
and these can be formulated as follows. If the restriction 
<f> c := 3?L_ C is nowhere-zero (as a function of x and y), 
then the normalized Higgs field $ c := $ C /|<I> C | is well- 
defined, and is a map from the torus T 2 to the 2-sphere 
S 2 . Here |$| 2 := — | tr ($ 2 ). So $ c has a degree (winding 
number) iV c G Z. If $ 2 is nowhere-zero for z > c, then 
by continuity the degree N z is independent of z, and is 
denoted N + ; similarly for N^. The boundary condition, 
motivated by the U(l) case, is: 

(c) D x $ -> and D y $ -> as z -> oo; if N + ^ 0, then 
|$|/jz — > 2-7T |JV+|, and \D Z $>\ is bounded, as z — > oo; 
if AT+ = 0, then |$| — > const and \D Z $>\ — > as 
z — > oo; and similarly as z — > — oo. 

We say that such a field has charge (— A^_, AT + ). 

For fields satisfying conditions (a)-(c), there is a topo- 
logical lower bound on the energy. The derivation of this 
is analogous to the one used for monopoles localized in 
R 3 . In order to get finite energy, we need to restrict to a 
finite cylinder —L<z<L] the condition (c) is adapted 
in the obvious way to become a condition at z = ±L. 
The energy density is 



£:=--tv[(D^) 2 + (B J ) 2 }. 



The first observation is that 



$5, 



dx dy, 



(6) 



(7) 



which is a standard calculation [6] . The energy is 

r-L 



E, 



-L 



dz I dx dy £ 



-~y tr(Dj$ + B 3 ) 2 d?x + J tY(Dj<S> ■ B 3 )d 3 x. 



This is obtained 0, Q by embedding the homogeneous 
Abelian field multiplied by a factor of -2, into SU(2): 



$ = 2mza 3 , A 3 = m(y, -x, 0)a 3 



(10) 



and then applying an SU(2) gauge transformation (which 
it is not hard to write down explicitly), in order to make 
the fields periodic. The energy density of (flTJ|) is easily 
read off, and is £ = 87r 2 ; so the Bogomolny bound ([8|) is 
saturated, as expected. 



IV. THE SU(2) CASE: NUMERICAL SOLUTION 

The problem is to investigate whether there is a SU(2) 
solution of the Bogomolny equations ^ with charge 
(0, 1) — and to see what it looks like. No existence 
theorem is currently available, and the best that can be 
done here is to use a numerical approach. The idea is 
to discretize the system as a lattice gauge theory, and 
to minimize the energy El numerically. A minimum 
which saturates the lower bound ([8|) , in other words with 
El — 8tt 2 L, should then be a solution. 

So we replace R 3 , or rather the region (0, 1] x (0, 1] x 
[-L, L], with a cubic MxMx (2LM +1) lattice T, so the 
lattice spacing is h = 1/M. The gauge potential is repre- 
sented, in the standard way, by assigning an element A of 
SU(2) to each link of T; the curvature is then an element 
B of SU(2) associated with each face. We represent the 
Higgs field by assigning an element <j> of SU(2) to each 
site of T; the covariant derivatives D<p are then elements 
of SU(2) associated with each link. The energy is given 
by the 'Wilson action' 

£J i = |^(l-trB) + |^(l-tr J D0), (11) 



faces 



links 



which is gauge-invariant on the lattice. 

The boundary conditions, which are also gauge- 
invariant, are as follows: 

• at z — L, we impose i tr<f> = cos(2-7r/iL); 

• at z = —L, we impose 6=1. 

These correspond, respectively, to the conditions |$| = 
2nL at z = L, and B = at z = —L, in the continuum 
version. 
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The initial configuration was constructed by starting 
with the lattice version of the homogeneous solution, 
namely 

<j> = cxp(2irihza 3 ) 
A x = exp(7ri%er 3 ) 
A y = exp(— 7ri/ixer 3 ) 
A z = 1; 

then gauge-transforming so as to make these periodic; 
and finally adjusting by setting <f) — 1 f° r z < 0, and 
interpolating so that A x and A y go to 1 as z goes from 
to —L. The resulting lattice field has charge (0, 1). 

Starting with this initial configuration, the energy (jTTJ) 
was minimized using a conjugate-gradient method, while 
maintaining the boundary conditions (and also, of course, 
the condition that A and <j> should be SU(2)-valued). 
This was done for various values of M (or equivalently 
h) and L. For M > 12, the results are already of high 
accuracy, in the sense that the value of El at its min- 
imum is within 0.5% of the Bogomolny bound, and it 
changes by less than that amount if M is increased. The 



Atx=1=y 
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FIG. 1: Charge, Higgs field and energy density of SU(2) 
monopole sheet 

graphs in Figure 1 depict the minimum-energy configu- 
ration with M = 14 (so the lattice spacing is h = 0.071), 
and L = 2. Its lattice energy is E L = 0.9997 x 8ir 2 L. 
The upper-left-hand graph shows the topological charge 
(winding number) N of $| , as a function of z; in par- 
ticular, we see that the solution does indeed have charge 
(0, 1). The norm of the Higgs field $ (at x = y = 1, 
and as a function of z) is shown in the upper-right-hand 
graph; as expected, it approaches a constant value as 
z — > — oo, as depends linearly on z as z — > oo. The lower- 
left-hand graph shows the energy density £ (or rather its 
lattice version) summed over x and y, as a function of z. 
We see that the energy density tends to zero as z — > — oo, 



approaches a constant value as z — > oo, and is peaked at 
z k 0.4. Finally, the lower-right-hand subfigure plots the 
energy density at z — 0.4, as a function of x and y. For 
\z\ > 1, the corresponding energy plots are essentially 
constant in x and y, but at the location of the monopole 
sheet (which here is at z ss 0.4) we see the individual 
monopoles (one in each fundamental cell of the lattice, 
of which four are shown) . 

V. CONCLUDING REMARKS AND OPEN 
QUESTIONS 

It is straightforward to change the x, y-periods by scal- 
ing: in fact, if {$, A,} is a solution with unit periods in 

x and y, then {$(r) := A$(r/A), A, (r) := AA,(r/A)} is 
a solution with periods A. More generally, it should be 
possible (although this has not been investigated) to con- 
struct analogous solutions in which the monopole sheet 
has a different lattice structure, for example hexagonal. 
(Compare the case of the Skyrme system, where a hexag- 
onal lattice is slightly more energetically favourable than 
a square lattice [l|, q}.) 

If one takes the gauge group to be SO (3) rather 
than SU(2), then the topological classification of doubly- 
periodic monopoles has an additional feature, since 
SO (3) bundles over the torus T 2 are not necessarily triv- 
ial — they are classified by Zg. Essentially, the conse- 
quence is that the charges (N-,N+) of SO(3) monopole 
sheets need not be integers, but may be half-integers. 
The simplest example Q is the homogeneous solution of 
charge (5,5), analogous to (fit))) , namely 

$<°> = mza 3 , 4° } = h*fa 0)^; ( 12 ) 
J 2 

which, in effect, lives on a non-trivial SO (3) bundle over 
T 2 x E. There should exist an SO (3) domain wall so- 
lution of charge (0, |) which is analogous to the charge 
(0, 1) solution. More generally, there should exist sheets 
of charge (p/2, q/2), where p,q € Z. If p and q are dis- 
tinct, then such a solution represents a wall between two 
distinct homogeneous phases; one of these could be the 
vacuum, if the corresponding charge is zero. 

The only 'visible' parameters in the numerical solu- 
tion presented above are those corresponding to transla- 
tions in x, y and z. (The position of the sheet, which 
turned out to be at z « 0.4, is in effect set by the bound- 
ary condition at z = L: in fact, by the fixed value of 
|$| at z = L.) Whether there are additional moduli is 
not known. For the homogeneous solutions, the moduli 
can be computed explicitly 0, Q : the solutions of charge 
(p/2, p/2), with p E Z, depend on 4p parameters; and the 
perturbations (normalizable zero-modes) can be written 
down explicitly in terms of theta-functions. (Only the 
p = 1 case was presented in [5] , but its generalization to 
p > 2 is straightforward.) This suggests that there might 
also be additional moduli in the general case of charge 
(pA?/2). 
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The main tool in the analysis of BPS monopoles has 
been the Nahm transform [H, [j| • The general pattern 
suggests that the Nahm transform of a monopole sheet (a 
solution on T 2 xK), will be another solution onT 2 xK- 
in other words, that monopole sheets are 'self-reciprocal' 
under the Nahm transform. The only current evidence 
in favour of this comes from the class of homogeneous 
solutions: the U(l) case was described in []|, and the 
SU(2) case is similar. 

Solutions (with prescribed singularities) of the Bogo- 
molny equations on C x /, where C is a Riemann surface 
and / is an interval, crop up in the area of supersym- 
metric gauge theory and branes (see, for example, [111]). 
The context of the present paper is rather different, but 
the solutions described here could be re-interpreted in 
D-brane language. 

There are many similarities between BPS mono- 
pies and Skyrmions (see, for example, [l2j ). In the 
Skyrme system, there is a wall-like solution analogous 
to a graphene sheet, and one may view fullerene-like 
Skyrmions as being shells constructed from such sheets 
Q; in fact, the sheet separates a region of 'Skyrmion core' 
from the region outside the Skyrmion. The monopole 
sheet described above appears at first sight to be rather 
different in nature, since neither side of the wall corre- 



sponds to the field outside a BPS monopole (which, of 
course, is that of a Dirac monopole). It may, however, 
turn out to be relevant as the wall of a magnetic bag [HI . 

To summarize: we have constructed, numerically, a 
solution of the Bogomolny equations representing a sheet 
of BPS monopoles: in general, it is a domain wall between 
two regions of constant energy density, either of which 
can be a vacuum region. 

Note added: I thank the referee for emphasizing that 
a purely-magnetic wall with vacuum on one side and a 
nonzero magnetic field on the other, cannot be static: it 
would accelerate in the direction of the magnetic field. In 
the BPS case, by contrast, the pressure on the wall from 
the Higgs field is exactly opposite to the pressure from 
the magnetic field, and so there is no obvious instability. 
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